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Some non-linear generalizations of classical Radon tomography were recently introduced by M. 
Asorey et al [Phys. Rev. A 77, 042115 (2008)], where the straight lines of the standard Radon 
map are replaced by quadratic curves (ellipses, hyperbolas, circles) or quadratic surfaces (ellipsoids, 
hyperboloids, spheres). We consider here the quantum version of this novel non-linear approach 
and obtain, by systematic use of the Weyl map, a tomographic encoding approach to quantum 
states. Non-linear quantum tomograms admit a simple formulation within the framework of the star- 
product quantization scheme and the reconstruction formulae of the density operators are explicitly 
given in a closed form, with an explicit construction of quantizers and dequantizers. The role of 
symmetry groups behind the generalized tomographic maps is analyzed in some detail. We also 
introduce new generalizations of the standard singular dequantizers of the symplectic tomographic 
schemes, where the Dirac delta-distributions of operator-valued arguments are replaced by smooth 
window functions, giving rise to the new concept of thick quantum tomography. Applications for 
quantum state measurements of photons and matter waves are discussed. 
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I. INTRODUCTION 

Quantum tomography has a long and interesting his- 
tory. The first intuition, hinting at a tomographic recon- 
struction of quantum states, dates back to 1933, when 
Pauli asked whether it is possible to uniquely associate 
quantum states with probability distributions, as in clas- 
sical statistical mechanics [il[2|- Pauli's observation was, 
in fact, more subtle and articulated, as he wondered 
whether two given position and momentum probability 
distribution functions are mathematically and physically 
compatible, and whether it is possible to unambiguously 
reconstruct the quantum state (the wave function) from 
their knowledge. The answer to Pauli's original question 
is negative 0, Ij] , but his general idea to associate quan- 
tum states with a set of probability distribution func- 
tions eventually led to the basic formulation of quantum 
tomography. 

The tomographic representation of quantum states 
is based on the Radon transform [al of the Wigncr 
quasiprobability distribution function Q. The applica- 
tion to the reconstruction of quantum states was pio- 
neered in the 80 's [7|, [SJ and led to the first experimental 
verifications in the 90's |9l-[ll|. Since then, the tomo- 
graphic approach to the analysis of quantum states has 
become a booming and consolidated field of investigation 
[la [jSj , leading to theoretical ideas and experimental 
proposals |14| - [23 |. thanks also to its important applica- 
tions in quantum information |23l . |24| | . 



The original Radon transform [5| maps functions de- 
fined on a two-dimensional plane onto functions defined 
on a two-dimensional cylinder. The key feature is that 
the transform is invcrtible. There exist several impor- 
tanl, generalizations of the Radon transform (sec, e.g., 
J25l - l27j ). More recent analysis have focused on symplec- 
tic transforms [23| , on the deep relationship with classical 
systems and classical dynamics 123, [SO] , on the formal- 
ism of star product quantization |3l| and on the study of 
marginals along curves that are not straight lines |32l | . 

The aim of this work is to study generalizations of the 
Radon transform to multidimensional phase spaces and 
to frameworks based on marginals along curves or sur- 
faces described by quadratic equations. We shall study 
applications to both classical and quantum systems. For 
classical systems we use the Radon transform of the prob- 
ability densities in phase space of a classical particle and 
build the corresponding tomographic map into Radon 
components of the initial probability densities. In quan- 
tum systems, we introduce a pair of operators, the den- 
sity operator and its Radon transform operator, which are 
related by an invertible map. In the context of Wigner 
function, these formulas provide an analogy of the Radon 
transform for classical probability distributions in phase 
space. 

The paper is organized as follows. In Section |ll] we 
review the symplectic tomographic approach to classi- 
cal mechanics. In Section IIIII wc provide a generaliza- 
tion of symplectic tomography that is based on marginals 



over ellipsoids and hypcrboloids and their boundary sur- 
faces, instead of planes and straight lines. The case of 
marginals based on circles and spheres with shifted cen- 
ters provides another type of interesting generalization. 
We also consider, by using a star-product framework, the 
generalization of Radon transform in operator form for 
quantum systems in Section ITVl Center of mass tomogra- 
phy [331 , symplectic tomography and its generalization to 
several quantum systems is discussed in Section |Vl using 
symplectic group Hamiltonians. A generic formulation 
of tomographic approaches that make use of Hamiltoni- 
ans which are linear forms in generalized Lie algebras 
is briefly summarized in this section. In Section IVII we 
consider thick symplectic quantum tomography and com- 
pare it with the thick Radon transform in Section IVIII 
Finally, some conclusions and perspectives are conveyed 
in Section EnH 



II. SYMPLECTIC QUANTUM TOMOGRAPHY 

Let /5 be a quantum state {p = p\p > Q). The quan- 
tum Radon transform (homodyne tomogram) is given by 



TZp{X,ip) = Tr/5 5 f XI— qcosip — psmip) , (1) 
with ip E [0, 2tt]. In terms of the Wigner function 

Wip,q)= J{q-^\p\q+^)d^, J W{p, q) dpdq = 27r, 

(2) 
the above expression reads 

TZpiX^ifi) = l.JdpdqW{p,q) 

xS {X — qcosip — ps'mip) . (3) 

The quantum symplectic [SJ] (or M^ [33 ) transform is a 
generalization of the quantum Radon transform ([1} and 
reads 



Wp{X,fi,iy) = (S{Xi-^iq-iyp) 



= Tr p S [ XI — ijq — lyp 



(4) 



where p and v are real parameters. The information con- 
tent of the two formulations is identical and is expressed 
by the relation 

yVp{X,rcosip,rsmip) ^ -npi—,ip] , (5) 

valid for any r > 0. Equation ^ is an easy consequence 
of the fact that the Dirac distribution is positive homo- 
geneous of degree —1. The two formulations ^ and ((4|) 
may differ in practice (the latter being easier to invert 
(35| ) and cease to be equivalent when the Dirac delta- 
function is replaced by a finite window function |36| . In 



the present article, we shall focus on the symplectic ver- 
sion ((4|). 

Equation Q can be rewritten as the "classical" tomo- 
gram of the Wigner function 



Wp{X, p,,y)^l- I W{p, q)5{X-pq- up) dpdq, 



(6) 



whose inverse transform reads 

W{j>, q)^^J Wp(^, f^, ^) e'(^-'''-''P) dX dp dy, (7) 
or, in terms of the density matrix, 

/5-7^ /Wp(X,/x,i/y(^^-^«-''P)dXd^di^. (8) 
27r J 

For completeness and for future convenience, we notice 
that the inversion formula of the Radon transform ([1} is 
easily obtained from Eqs. (O and ((8|) and reads 

-I f /"OO flTT 

p = ^J ^^ J dry d^TZpiX,^) 



^ir{XI-~-q cos ip~psin ip) 



(9) 



The generalization of ([6]) to the multimode case can 
be achieved in several ways. Let us consider the case of 
two modes for simplicity. One can either introduce the 
ordinary tomogram 



^piX,p,iy) 

S{Xii - piqi - vipi) 5{X2i - ^292 - V2P2) 



(10) 



where X = (Xi,X2), p = {pi,p2) and u = (1^1,1^2), or, 
alternatively, the "center-of-mass" tomogram 



WcmiX,p,iy)^(S{Xi- piqi - vipi - P2q2 - i^2P2) 



(11) 

The Weyl map ^ provides again the possibility of ex- 
pressing pUj) and pT|) as "classical" tomograms 



Wp(X,A^,z7) 



1 



dpdqW{p,q) 



5 {Xi - piqi - vipi) S {X2 - P2q2 - ^^2^2) , 

(12) 



WcmiX, p,lJ) = 



1 



dpdqW{p, q) 



(27r)2 

5{X - piqi - Vipi - p2q2 - V2P2) 



(13) 



where W{p, q) is the Wigner function of the two- mode 
system, normalized to (27r)^. 



The inverse transform reads 



Wip, q) = / Wcm(X, M, y) SX^fi-^-u-p) 



dX dfl dv 



for the center-of-mass tomogram, and 



(14) 



(15) 
for the ordinary symplectic tomogram. The inversion 
formulae for the density matrix are 

p = J W,iX, p, V) ci(x.i+x.i-,.,^-..^) d^^^ (Ig) 
for the symplectic case, and 

P = / Wc„.(X, M, V) <}^-^-^-'-^ ^^ (17) 

for the center-of-mass tomography. The advantage of the 
symplectic tomogram is that it permits the splitting into 
subsystems for separable states, which is not possible in 
the center-of-mass tomogram. This property might be 
very relevant for the fast reconstruction of non-entangled 
states. 

In the case of M subsystems with their own centers of 
mass, tomography is defined by 

WcmiXi, pi,Vi,X2, P2, 1^2, ■■■ ,Xm, f^M,VM) 
I M \ 



\{5[x,l-p,-l-iJ,-T>,)),{l%) 



where fli = ip-ii,Pi2, ■ ■ ■ ^fJ-iN,), i'i = {^ii^i^i2i- ■ ■ i^iN,), 
with i = 1, 2, . . . , M and 



M 



Y.N,=N. 



(19) 



M X M symmetric operator B and two M-dimensional 
vectors x and a, let 



X ~ {x — a) ■ B{x — a). 



(20) 



where x ■ y denotes the scalar product of vectors x and 
y. The classical tomogram of a function f{x) reads 

ujf{X,d)^ / f{x)S{X ~{x~d)-B{x-d))dx. 



(21) 

Observe that the ^-function is supported on the quadrics 
(ellipsoids, hyperboloids, etc., that can also have degen- 
eracies) defined by Eq. (j20p . The inverse map is 

fix) = ^H^ / dXddujfiX,d) 



xe 



i{X-{x-a)-B{x-a)) 



(22) 



The solution of the quantum problem is similar. With 
a slight abuse of terminology, we shall speak of multidi- 
mensional phase spaces and quadratic "Hamiltonians" of 
the type 



H=^Q-BQ + C- 



(23) 



where 



Qj = Pj, Qn+j ^qj, (j = 1, . . . , iV), (24) 

are the momentum and position operators, and C and B 
denote an 2Af-dimensional vector and a 2N x 2N sym- 
metric matrix, which parameterize the different types of 
the quadratic Hamiltonian. Observe that since we are 
working in a symplectic framework, M — > 27V. 

The quantum counterpart of the classical tomogram 
(PT|) reads 



The inversion formula is given by 

P = yVcm{Xi,fliyi,X2,P2,'^2, ■ • ■ , Xm, pm, vm) 

J E, {xii-P, ■?, -P. -P- ) dXdpdiJ 
(27r)*^ ■ 

The symplectic tomogram can also be constructed in the 
case of multimodal systems. 

III. QUADRATIC TOMOGRAMS 

One of the main objectives of this article is to 
study generalizations of the quantum Radon trans- 
form to marginals along curves or surfaces described by 
quadratic, rather than linear equations. The solution of 
the classical problem was obtained in Ref. [S^] ■ Given an 



WiX,p,u) 



6 XI 



•^nij 



(25) 



where 



^fiu = :^{Q-r)-B{Q-r) + C-{Q-r) (26) 

and f = [p, V) is a 2A^-dimensional vector. The coordi- 
nate X plays the role of "energy" . The Hamiltonian can 
be degenerate, depending on the type of eigenvalues of 
the 2N X 2 A^ symmetric matrix B [3a] . The parameters p 
and V have the meaning of shift parameters of the centers 
of the operatorial quadratic curves (surfaces). 

From the Hamiltonian (|26p one can derive a quantum 
tomographic map similar to that obtained in the classical 
case, by using the Wigner function: 



y^{x,p,v) 






{fi,)W{q,p), (27) 



where the Hamihonian Mpp is given by ([25]) with the op- 
erators q and p replaced by c- numbers q and p. Since Eq. 
^7\ is equivalent to Eq. ((2T|) . the reconstruction func- 
tion follows directly from Eq. (|22|) with the replacement 
B — > ii3, namely, 

P = />V(X,^,.-)eH^S-%^^)|detS|^i^^. (28) 

The generalization to multipartite systems is straight- 
forward. If we have M subsystems we can define the 
multipartite tomogram 



M 



Wem(^, P, V)= W{q,p) Yl 6 (X,i - %^,- 



J = l 



^Pjdqj 

(2^) 



Ni 



X = {Xi,X2,...,Xm), P^ifij). v = {v,). (29) 
whose inverse transform is 

» M 

p = /we.„(x,M,^)n^'*''^'"*"^'^-^ 



X |detB 



dXj ApLj Avj 



(30) 



IV. STAR-PRODUCT FRAMEWORK 

It is interesting to see how the transform with 
quadratic curves in phase space can be written in the 
form of a star-product quantization [3l|. Any star- 
product scheme is aimed at constructing a bijective map 
of operators A acting on a Hilbert space % and the space 
of functions /a defined on a manifold A^ . The map can 
be constructed in terms of two operator- valued functions 
D and U defined on M and called "quantizers" and "de- 
quantizers" , respectively. These two functions satisfy the 
identity 



Tv{b{x)U{y)] ^6{x-y) 



(31) 



for every x,y G A4. The map is then defined by the 
formula 



(32) 



(33) 



The star product is defined on the space of functions /a 
by 



and the inverse map is defined by 

A^ j fA{x)b{x)dx. 



fA^fs — /aB, 



(34) 



The kernel of this associative star product is then given 
by 

K{xuX2,X3) ^Tr (^b{x,)D{x2)Uix3)'^ (35) 



and satisfies 

fA-*:fB{x)= fA{xi) fA{x2)K{xi,X2,x) dxidx2- (36) 

The key observation is that the tomogram can be in- 
terpreted in terms of a dequantizer 

(37) 



u{x,^i,,;)^s{xi-mpp 

and a quantizer 



D{X,fl,i7) 



\dctB\ 



,i(xi- 



(38) 



Thus, for an arbitrary observable A, one can introduce 
the quadratic tomogram 



(39) 



WAiX,fl,v)^TiAS{XI- 



'^liv ) 1 



Also one can reconstruct the observable operator from its 
quadratic tomogram 



A 



Wa{X, fl, v) b{X, fl, V) dX dildv. (40) 



There exist also tomograms based on shifts of quadratic 
curves in phase space, but, as in Eq. Q for symplectic 
tomograms, we used instead rescaled position and mo- 
mentum operators. 

One can combine both types of tomograms, derived by 
shifts and rescaling of variables, obtaining a new class of 
tomograms, which are the quantum version of the classi- 
cal tomograms defined by j32l | 



W/(X,^,l7): 



S{X - fl-q- v{q,p)) f{q,p)dqdp, 



where p and q are vectors in R" and 

n 



(41) 



(42) 



This map corresponds to a deformation of the standard 
multidimensional Radon transform by means of the fol- 
lowing diffeomorphism of R^"-\ [^Aiq^p) ■ qj = 0} 

ilt^Pj)^ {xt , yj ) = iq^ , qjPj ) , (43) 

whose Jacobian is 

d{x, y) 



Jiq,p) 



d{q,p) 






(44) 



The inverse map is given by 
f{q,p) = 



dX dfl dv , , ^ ^s 






i(X-p,-q-v{q,p)) 



(45) 



This corresponds to the highcr-dimensionalgeneraUza- 
tion of the Bertrand-Bertrand tomography [3] . 

The quantum extension is straightforward. One intro- 
duces the Hamiltonian 

N N 

^ii^{q,f) = Yl ^'^' + 2 ^ ^j(9iW +Pj^j)- (46) 

The dequantizer operator 

U{X, C V) = 5(XI- %(g-, p)) (47) 

yields the quantum tomogram 

W(X, C, y)^i:ip5 [xi - i^^Mp)) , (48) 
and corresponds to the quantizer 



N 






\{Xl-n.^(q,p)) 



(49) 



This permits to recover the quantum state from its to- 
mogram. The modulus of the operator is defined by the 
Weyl quantization of the symbol \qj\. 



V. GROUP HAMILTONIANS 

The applications considered so far involve position and 
momentum coordinates. In the linear and quadratic case, 
the Hamiltonians from which we constructed the quan- 
tizers and dequantizers can be written in the form 



Hg^ — 7 \go)aLa, 



(50) 



where La arc generators of Lie groups. In the examples 
considered in this article, they belong to representations 
of the inhomogeneous Lie group ISp(2A^, M) and its sub- 
groups. The group ISp(2A^, M) is defined as the group 
of affine transformations of M^^ = T*W^ which preserve 
its natural symplectic form. It the semi-direct product 
of the translation group and the linear symplectic group. 
The reason for this choice is because this group acts on 
the Heisenbcrg-Wcyl group as a group of automorphisms. 
Within this framework, the delta function is expressed as 

5 (xi-Y,{go)aL^ ^ i-y'die'*(^i-5:»(^o)"^") 

(51) 
and 

;7(.g(i)) =e-'*(^»(9°)°-^°) (52) 

is the operator representation of the Lie group. 



Thus, the tomogram can be rewritten in a form that 
only depends on group element parameters [go) a and the 
Radon variable X: 



W(X, (go)a) 



5[xi-Y,{9o)cX 



Tr / dtpc''^U{g{t)) 



die'*^Tr 



\(pU{g{t))). (53) 



Therefore, the tomogram is related to the orbit of the 
group and the Fourier integral of the trace of the orbit 
in the group representation. This group-theoretical rep- 
resentation permits the extension to more general cases 



VI. THICK QUANTUM TOMOGRAPHY 

We now turn our attention to thick tomographic maps 
j3a | , which is a more realistic approach for practical appli- 
cations, because instead of marginals defined over lines, 
as in the classical Radon transform |5| or quadrics as 
in the quadratic generalized Radon transform [30|, l32l |. 
it involves a thick window function S. This is convo- 
luted with the tomographic map and concentrates the 
marginals around some given background curves (that 
can be lines or quadrics), without resorting to a singu- 
lar delta function. For example, if the weight function S 
is a step function, it defines marginals along thick lines 
or thick quadratic curves. In the quantum case, this 
amounts to replacing in the definition of the dequantizer 
U{x) the Dirac delta-function by the weight function S, 



U{X, [g 



°)„)-sfxi-^(.g„)„i„j 



(54) 



For the symplectic quantum tomography, one has the 
dequantizer 



(55) 



U{X, n,!^) = E{XI- fiq-iyp 
The new tomogram reads 

WsiX, ^i, iy)=TTpE{XI- pq - vp 



(56) 



Using the Weyl map one obtains a tomogram for the 
Wigner function 

WsiX,p,iy) = j- J Wip,q)EiX - fiq-iyp) dpdq. 

(57) 
The interesting property of the above formula ((57)) is that 
it can be inverted in complete an alog y with the classical 
thick tomography introduced in [36[. The thick tomo- 
gram can be expressed in terms of standard symplectic 
tomograms via a convolution formula 



W. 



:{X, p,iy) = J W(y, p, iy)E{X- Y) dY, 



(58) 



which leads to the exphcit construction of the inverse 
transform. Indeed, the inverse transform is obtained by 
means of a Fourier transform of the convolution integral 



W{p,q) - ^ / WsiX,fi,iy)e'^^-'"'-'P^ dXdfidiy, 

ZTT 



where 



Ms = 



1 



(59) 



S(-l)= / S(z)e'"dz. (60) 



S(-l) 
In invariant form, the state reconstruction is achieved by 

p^-^ fw^iX,n,iy)c''^^^-'"i~''P^ dXd^di^. 
27r J 

The quantizer operator in thick symplectic tomography 
is 



2tt 



(61) 



Let us now consider a particular example of thick tomo- 
gram to illustrate the potentialities of the new method. 
If the weight function is a gaussian function 



S(z) 



1 



V27rcr2 



(62) 



which tends to the delta distribution in the ct — > limit, 

(63) 

the thick tomogram of the coherent states |Q;)(a| read 



lim^(2:) =(5(z), 

CT— S-O 



W,"(X,/x,^) = ^===i===e Ti^ff^, (64) 

where 

X = y2^Rea + \/2i/Ima. (65) 

For the vacuum state |0)(0|, the tomogram reads 

1 



wr{x,p,,,) = 



The quantizer reads 



e ^'"+-"+2-^ . (66) 



ATT 

and the dequantizer is given by 



Ua{X,p,u) 



1 (Xl-l.Lq-up)'^ 



V27rcr2 



(67) 



(68) 



One interesting property, that is preserved by the 
smoothing of the tomogram, is that the marginals 
We.{X,Ii,v) are also probability distributions. In the 
limit cr ^ 0, S(z) -> 5{z), S(-l) = 1, TVs = 1. 



In the case of multimode systems with quadratic 
Hamiltonians, the dequantizer and quantizer of thick to- 
mography with a Gaussian weight function S become 

U,{X, fl,,7) = E (xi - 1(0 - f) • S((^ - f)) (69) 

and 

D4X,fl,^) = IdctBl i(A-i-|(Q--.).B(9--o) 

^ '' ' (27r)^+iS(-l) 

(70) 
respectively. Thus, the quantum thick tomogram of a 
multimode state with weight function S reads 

Wh(^,/I,?) = [ ^W{q,p) 



(27r) 



1 



El^X--{Q-r)B{Q-f)j, (71) 

in terms of the Wigncr functions of the state, and the 
inverse formula 



(27r)^+iS(-l) 

gi(jf-i(Q-r)B(Q-r)) 



(72) 



permits to reconstruct the Wigner function from its to- 
mograms. 

The proof of the inverse formula follows the same steps 
of the case of linear symplectic tomograms. Indeed, if we 
plug ([7T|) into the right hand side of ([7^ we get 



|detB| /■ ^^ j^^_^ f dp' dq' 



(27r)^+iS(-l) 
|detB| 



dX dfl dv 



(27r 



,N 



W{q\p') 



^(x- i(Q' - r^B{Q' - r)\ ci^-k(Q-rlB{Q-r)) 

dXdYdfldu [ ^^^Wiq",p') 



(27r)^+iS(-l) 



(27r) 



S{Y ~^iQ'- r)B{Q' f))E {{X Y)l) 

gj(y-i(Q-r)i3(Q-7^) ^iiX~Y) 

\dctB\ /■ ,7 ,^, ,^/„.^^, -.,. f dfldv _ ^z 

^-^^J dZdpdqWiq,p)J J^-'-^" 
^i{rB(Q-Q')-k_QBQ + ^_Q' BQ') 



■S(Z)e' 



J— W{q,p) I f- S(Z) e'^ = W{q,p) 
S(-l) J 27r 



(73) 



Let us consider the example of a quadratic thick tomo- 
gram governed by a Hamiltonian with 



B = 



l2,0 
0,l2 



(74) 



and Gaussian weight function So- for a system with one 
degree of freedom. In that case, one might think that the 
tomogram is centered on a quadratic curve 



X = \{p-l^f + \{q-y)^. 



(75) 



The tomogram is given by 



VIII. CONCLUSIONS 



Waix,^l,l^) 



dpdq W{q,p) 



27rV27rcr2 



(76) 



in terms of the Wigner function and the inverse transform 



e 2 



w{q,p) = -— / dxd^ldlyw^{x,^i,l^) 



(77) 



allows a reconstruetion of the Wigner function from its 
thick quadratic tomograms. 



VII. COMPARISON WITH THE THICK 
RADON TRANSFORM 

A comparison with the quantum Radon transform ([1]) 
is useful and clarifies the merits of the approach taken in 
this article. Let us consider the thick version of ([T]) , with 
a window function S 

nE{X,ip) =Trp E{Xl-qcosip-psm^). (78) 

Now, in general, relation ([5]) (a consequence of the pos- 
itive homogeneity of degree —1 of the Dirac distribu- 
tion) does not hold anymore for TZb. and Wh- As a con- 
sequence, one has to deconvolve 7?.h — whenever this is 
possible — in order to get TZp and thus p. Since 



7^H {x,ip)^ jE{x-Y)np (y, ^p) dy. 



(79) 



one gets 



dX7^H(X,v3)e-"=-^ =S(fc) / dXnp{X,ip)c-'' 



ikX 



(80) 
Equation (|80p can be inverted for any TZp if and only if 

S(fc) ^ for all k. In this case, from ([50)1 and © one 
gets 



1 



27r 



/5 - 7^ / dX / dr / d^ 

^0 

ir{XI—q cos ip— psin (f) 



OO n27T 



7^5(x,^) 

S(-r) 



xre 



(81) 



Therefore, in contrast with the thick symplectic trans- 
form, the thick Radon transform cannot be inverted 
for arbitrary window functions. Moreover, even when 
the inversion is possible, one needs to have a complete 
knowledge of the window function in order to deconvolve 
the thick homodyne tomogram and perform its inversion 

m. 



There are many possible interesting applications of the 
quantum tomograms introduced in this article by using 
quadrics and window functions in the phase-space. The 
measurements of variables that are quadratic forms of 
positions and momenta are equivalent to measuring the 
energy of vibrations with changing parameters, like min- 
ima of potential energy, elastic forces, and so on. This 
procedure can be used to measure the state p of trapped 
ions by varying the parameters of the trapping potential. 

Homodyne detection, by using non-linear crystals as 
beam splitters, also creates quadratic forms of the photon 
quadrature components which, in principle, can be mea- 
sured by varying the parameters of the quadratic forms. 
Thus, tomography with the generalized quantized ver- 
sion of the Radon transform provides new possibilities to 
measure quantum states. 

From a formal viewpoint, the new tomographic maps 
and their inverse can be formulated within the frame- 
work of the star-product scheme, with quantizers and 
dequantizers for multidimensional systems. Some of the 
constructions involve generic quadratic forms of the dif- 
ferent position and momentum operators. The inhomo- 
geneous symplectic group of the corresponding multidi- 
mensional phase spaces and its subgroups play a signifi- 
cant role in the formalism. The new quantum tomograms 
can be considered as the quantum version of the general- 
ized classical tomograms introduced in our previous work 
[32| and involve quadratic surfaces on phase spaces. The 
correspondence between both approaches is established 
by means of the Wcyl map. 

In quantum optics, the new non-linear Radon trans- 
forms can be easily extended to the quantum domain 
by using the Weyl- Wigner map. The results of this ar- 
ticle show that the reconstruction of the Wigner func- 
tion using optical or symplectic tomography based on 
straight-line Radon transform can be extended to situa- 
tions in which the marginals in phase space are measured 
for curved hyperbolas or ellipses. In particular, parabolic 
tomography could be implemented with the recently ob- 
served accelerated Airy beams 



Finally, the other generalization introduced in this pa- 
per, named "thick" quantum tomography, entails that 
the effect of a non-singular window function not only 
smooths the tomograms, but in fact matches better all 
practical applications, where thick tomograms arise in a 
natural way. The striking fact we pointed out in this work 
is that the reconstruction procedure in this case is almost 
independent on the window functions, with an enormous 
advantage due to the fact for most experimental devices 
the window function is not completely known. 
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